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Abstract 

We present a numerical method for investigating the non-perturbative quan- 
tum mechanical interaction of light with atoms in two dimensions, without 
a basis expansion. This enables us to investigate intense laser-atom inter- 
actions with light of arbitrary polarization without approximation. Results 
are presented for the dependence of ionization and high harmonic generation 
on ellipticity seen in recent experiments. Strong evidence of stabilization in 

circular polarization is found. 
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In the last ten years the interaction of intense laser light with atoms has been of great 
interest because of the many novel non-perturbative phenomena that have been observed, 
such as Above Threshold Ionization (ATI), High Harmonic Generation (HHG) and Stabi- 
lization Jlfl. The strongly non-perturbative nature of all these phenomena makes theoretical 
analysis difficult. However, considerable success has been made in using semi-analytic ap- 
proaches and numerical methods for integrating the Schrodinger equation. The numerical 
approach has been used extensively in the past ten years due in part to the increase in 
available computing power, the lack of approximations (except for those inherent to the 
Schrodinger equation and the dipole approximation) and the rather limited applicability of 
analytic approaches. Despite the success of the numerical methods most investigations have 
dealt with only linearly polarized light: the reason for this is that at least two dimensions are 
needed in the case of arbitrary polarization, so established one-dimensional approaches || 
are inappropriate. Also, current three-dimensional approaches are unable to deal with 
the enormous basis set needed to represent adequately the evolving wave function unless 
a parallel supercomputer [7| is used. The approach used by many authors || to describe 
atomic dynamics in strong laser fields expands the atomic wave function into sums of prod- 
ucts of radial wave functions (channels) and angular factors. These channels are coupled 
together by the very strong laser field, and for linear polarization, the channel couplings 
lead to tri-diagonal sparse matrices amenable to efficient numerical solution. This approach 
entirely fails for polarizations other than linear, as the matrix sparseness properties totally 
disappear ||. 

Recent experiments on the very strong dependence of high harmonic generation 

on polarization ellipticity demand theoretical approaches to intense field ionization going 
beyond the existing methods which rely on linear polarizations. In this letter, we introduce 
a method for integrating numerically the two dimensional Schrodinger equation without 
using a restrictive basis expansion. This enables us to study light of arbitrary polarization 
and realistic frequencies, interacting with an atom. Fully quantum mechanical results for 
ionization and high harmonic generation are presented for circular and elliptically polarized 
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light, showing the strong ellipticity dependence seen in recent experiments @,[T(J. Also, 
we show that in the intermediate regime between tunneling and multiphoton ionization the 
ionization dependence on ellipticity is opposite to that seen with purely tunneling dynamics; 
this may be critical for experiments towards generating attosecond pulses. 

As in all studies of this kind the starting point is the Schrodinger equation (atomic units 
are used throughout): 
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where ^>(x,y,t) is the two dimensional wavefunction, /(£) the pulse envelope and e the 
ellipticity. The second term on the right hand side is the smoothed Coulomb potential 
used to avoid numerical problems associated with the singularity at x = y = (a = 0.8 so 
that the ground state has the same binding energy as hydrogen, i.e. —0.5 a.u.). The pulse 
envelope, f(t), used here is a 2-cycle linear turn-on followed by 2 cycles of constant intensity 
and then a 2-cycle linear turn-off, which ensures no residual drift motion in the free electron 
displacement. The ellipticity is defined such that e = represents linear polarization along 
the x-axis and e = 1 circular. 

Equation [l| is integrated using a method analogous to that implemented by Grobe and 
Eberly |T]J and others in their studies of one-dimensional two-electron systems in which 
both electrons are restricted to one dimension. In our case there is only one electron with 
two degrees of freedom. The core of the approach is the split step method [12|] in which the 
time evolution operator is split into kinetic and potential operators: 

U(At) = exp[-iTAt/2] exp[-iVAt] exp[-iTAt/2] + 0(At 3 ), (2) 

where T is the kinetic part of the Hamiltonian (including the dipole interaction, in a 'minimal 
coupling' form) and V the potential part. The action of the kinetic operator is efficiently 
carried out in Fourier space, while the action of the potential operator is carried out in real 
space. An absorber is used to remove any part of the wavepacket reaching the boundaries, 



so that artificial reflections are avoided. It is assumed that any wave function reaching the 
end of the grid box represents ionization, and so the remaining norm gives the necessary 
information about the amount of ionization occurring. This means that any population in 
high lying Rydberg states may also be removed; however, this is negligible for the parameter 
regime of interest here, as we have verified by varying the box dimensions. In all the cases 
presented, the initial wave function is the ground state, found by imaginary time integration. 



A more detailed discussion on various computational aspects can be found in ref. [JTT . 

In Fig. 1 we show the wavepackets generated after 3.5 cycles of evolution for the previ- 
ously described pulse shape with laser light of wavelength 526 nm (u = 0.0867 a.u.), i.e., the 
wavelength of a frequency doubled Nd:YAG laser and intensity 3.51 x 10 lA Wcm~ 2 (0.01a.u.). 
For the case of linear polarization it is clear that the wavepacket is concentrated about the 
polarization axis, but there is a non-negligible width due to transverse wavepacket spreading. 
Interestingly, side lobes can be seen in the wavepacket; these are formed by a combination 
of the slow transverse spreading, the relatively fast laser driven motion in the polarization 
direction and the Coulomb attraction (rescattering) . 

The behaviour of the wavepacket created by circular polarization is very different. The 
packet resembles a spiral with modulations within the tail, formed as the wavepacket streams 
away from the atom down the potential barrier in the radial direction. As well as being in- 
teresting in their own right the wavepackets show pictorially that the recollisions necessary 
for efficient HHG can only occur for light with e ~ 0. The complex structure of the wave 
function in the linear polarization case reflects the fact that rescattering from the atomic 



core, according to the recollision model of HHG [|T3| p~5|| , produces strong time-dependent 



interference patterns in the vicinity of the atomic core, which are responsible for much of 
the structure in the spectrum of HHG Those interferences are absent in the circular 
polarization case, as is clear from Fig. 1(b), due to the fact that the probability of rescat- 
tering with the atomic core is negligible. We have observed a smooth transition from the 
recollision behaviour resulting from linear polarization to the swirling tail of the circular 
polarization, by integrating the wavepacket dynamics for ellipticities ranging from to 1. 



This fact can be shown explicitly by investigating the dependence of ionization and HHG 
spectra on ellipticity. This is shown in Fig. 2: in (a), for a fixed intensity of 0. 01a.it. the 
normalization at the end of the pulse is plotted as a function of ellipticity, where the results 
have been scaled such that the normalization for the linear case (e = 0) equals unity. One 
has to keep in mind that since the intensity is kept constant for varying ellipticity, the am- 
plitude of the electric field in the initial (linear) polarization direction changes. The reduced 
ionization in the initial direction is compensated by the possibility of electron ejection in 
all other directions. The result is a net increase of total ionization yield, especially for such 
high intensities that Over the Barrier Ionization (OBI) occurs. This effect is even greater if 
equal electric field amplitudes are compared rather than equal intensities. In this case the 
probability of ionization is obviously enhanced, as the electron can escape from the atom 
by tunneling through or over the barrier at any time in the direction of the instantaneous 
electric field. 

Clearly, for these parameters we find that as the ellipticity is increased the amount of 
ionization increases by a factor of approximately 3. Past experiments [171,18 have found that 



the threshold for ionization is higher for circular polarization. However, these experiments 
have been carried out in the tunneling regime where standard tunneling theory [19] predicts 



higher ionization rates for linear polarization. The results presented here are not strictly 
within the tunneling regime (due mainly to the relatively high frequency) and represent new 
behaviour. In Fig. 2 (b), we have used the same parameters as in Fig. 1 (the spectra are 
rescaled such that the fundamental peak is equal to one in intensity for each ellipticity), and 
now find that the harmonics act as seen in experiments PJT0"||: as the ellipticity is increased, 



the harmonics become less intense until no harmonics are discernible from noise. Only the 
first 4 harmonics can be seen for e = 0.6, and the peaks are found to be rather broad, 
indicating that they are generated on a short timescale. In addition there seems to be a 
correlation between the ionization behaviour and the ellipticity beyond which the harmonics 
are destroyed: the normalization flattens out after approximately e = 0.7 which is also the 
ellipticity at which the harmonics disappear. 



To highlight the very different behaviour of the ionization for various ellipticities we have 
calculated the normalization at the end of the pulse for linear and circular polarization as 
a function of intensity (Fig. 3). At low intensities we see the expected perturbative multi- 
photon behaviour, observed in experiment ||10|| , that linear polarization produces a greater 
amount of ionization because of the transition rule restrictions on circular polarization tran- 
sitions. However, as the intensity is increased to the critical intensity for over the barrier 
ionization (around 0.005a.w. of intensity) we observe that the behaviour changes such that 
circular polarization becomes more efficient in producing ionization and indeed for much 
of the intensity range shown it is at least an order of magnitude greater. Note that the 
normalization curve for linear polarization is not smooth: this was previously observed by 
Pindzola and Dorr |2(| and LaGattuta [^1|] and assigned to AC Stark shifted resonances 



and threshold effects. However,we find, as in ref. |2T| that for circular polarization these 
structures do not occur, perhaps because of the more restrictive circular dipole selection 
rules, at least in the perturbative regime of intensities. 

Finally, we observe a minimum in the normalization implying that stabilization is oc- 
curring. This is in agreement with calculations made by Pont and Gavrila JZIJ for circular 
polarization using time-independent high frequency Floquet theory. The lowest frequency 
which they have used is u = 0.125 a.u., while we have found stabilization, albeit not to the 
same degree, for a frequency 1.5 times smaller. Interestingly, for the same parameters linear 
polarization shows no signs of stabilization. This may be due to the atom not surviving 



the passage through 'death valley' f23| . We are currently investigating adiabatic stabiliza- 



tion with much higher frequencies and intensities, the results of which will be presented 
elsewhere. 

In summary, we have investigated numerically, with two-dimensional ab initio calcula- 
tions, the dependence of ionization and HHG on ellipticity. In particular, we have found 
that the harmonics are destroyed by elliptical light, as seen in recent experiments [p|,[10|] . 
In a parameter regime where the tunneling picture is not entirely valid, circular polariza- 
tion produces a greater amount of ionization than linear. We have also found evidence 

6 



for stabilization in circularly polarized light in qualitative agreement with previous time- 
independent calculations. The method we propose opens up the possibility of studying fully 
quantum-mechanically wavepacket evolution in time-dependent ellipticities proposed for the 



generation of ultrashort pulses of harmonics P^| . As well as this exciting prospect, the 
unified description of the strong-field phenomena, based on the exact wavepacket dynamics 
presented here, will open the way to resolve problems associated with, e.g., the accuracy 
of the 'two-step' model for harmonic generation, time- dependent laser polarization, con- 
trol of harmonic generation in time, ionization with elliptical light, elliptical stabilization of 
wavepackets, angular distribution of photoelectrons in ATI, and many other important ques- 
tions in strong-field physics that could be addressed only by approximate and speculative 
theoretical models in the past. 
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work has been funded in part by the UK Engineering and Physical Sciences Research Council 
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FIGURES 

FIG. 1. Snapshots of the probability density taken after 3.5 cycles for (a) linear polarization 
and (b) circular polarization with intensity 3.51 x W 14 Wcm~ 2 and wavelength 526 nm. Distances 
(X,Y) are shown in atomic units (a.u.). 

FIG. 2. (a) Rescaled normalization at the end of the pulse as a function of ellipticity for a 
6 cycle pulse with maximum intensity 3.51 x W w Wcm~ 2 and wavelength 526 nm. (b) Rescaled 
HHG spectra as a function of ellipticity for the same laser parameters. 

FIG. 3. Normalization as a function of peak laser intensity for a 6 cycle pulse with wavelength 
526 nm, for linear and circular polarization. 
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This figure "figla.gif" is available in "gif" format from: 



http://arXiv.org/ps/quant-ph/9708033vl 



This figure "figlb.gif" is available in "gif" format from: 



http://arXiv.org/ps/quant-ph/9708033vl 
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